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Abstract. We explore planet formation in the Quark Nova scenario. If a millisecond pulsar explodes as a Quark 
Nova, a protoplanetary disk can be formed out of the metal rich fall-back material. The propeller mechanism 
transfers angular momentum from the born quark star to the disk that will go through viscous evolution with 
later plausible grain condensation and planet formation. As a result, earth-size planets on circular orbits may 
form within short radii from the central quark star. The planets in the PSR1257+12 system can be explained by 
our model if the Quark Nova compact remnant is born with a period of ~ 0.5 ms following the explosion. We 
suggest that a good portion of the Quark Nova remnants may harbour planetary systems. 
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1. Introduction 

The first planetary system beyond our own was found or- 
biting the millisecond pulsar PSR 1257+12 (Wolszczan 
& Frail 1992; Wolszczan 1994). There are observations of 
three planets, Mi = 0.015M ffi / sinii, M 2 = 3.4M®/ sinii 
and M3 = 2.8M0/sini3, with the semimajor axes r\ = 
0.19 AU, r 2 = 0.36 AU and r 3 = 0.47 AU, respectively. 
Here is the inclination of the orbit with respect to the 
observer. The planets have very small eccentricities, i.e. 
their orbits are practically circular. The nature and ori- 
gin of these planets is still unknown and debated (see e.g. 
Lin et al. 1991; Phinney & Hansen 1993; Podsiadlowski 
1993; Miller & Hamilton 2001). Models can be divided 
roughly into two classes depending on whether the planets 
are formed before or after the supernova (SN) explosion. 
In the pre-SN scenarios the planets form in a planetary 
disk before their host star exploded as a supernova. In the 
PSR 1257+12 system the three innermost planets would 
have been inside the envelope of the progenitor star be- 
fore the explosion, evaporating or spiraling to the center of 
the star. Planets surviving the supernova explosion could 
have been far away from the star, and they could have 
migrated to their present distances. It is argued that this 
scenario is not very probable (Miller & Hamilton 2001). 
It has also been argued that the explosion would unbind 
planets initially present, and any remaining planets would 
have high eccentricities (Podsiadlowski 1993). 

More attention has been focused on post-SN scenarios. 
Summarizing the work of Miller & Hamilton (2001), it was 
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found that if the planets form before or during the possible 
spin-up period of the pulsar (following accretion from a 
companion star which may also provide the disk material) 
they will be evaporated by the accretion luminosity unless 
the disk is thick enough. If the planets are formed after 
the spin-up they must be created from some remnant disk, 
but then the particle luminosity from the neutron star 
(NS) is sufficient to disperse a tenuous disk material faster 
than it can be supplied. Added to that, the disk must be 
optically thick to the high energy particles of the pulsar 
wind as to allow for planet formation. That is, the disk 
must be massive enough ~ 10 28 g at a radius of 0.19 AU 
(the radius of the innermost planet in the PSR 1257+12 
system; see Miller & Hamilton 2001 for more details). The 
SN recoil scenario (Phinney & Hansen 1993) where the 
NS receives a kick in the direction of the companion as 
to capture enough matter from the companion is favored. 
This scenario also implies that many isolated NSs (maybe 
up to around 10 %) may have planets in orbit around them 
and must thus await further observational evidence. 

In general, it seems that very special kind of circum- 
stances are required to produce planets around a pulsar. 
Any model must demonstrate that the disk can expand to 
the appropriate radius with suitable conditions for planet 
formation on reasonable time scales (Bodcnheimcr 1993). 
Furthermore, the mechanism cannot be too efficient, or 
otherwise more planetary systems would have been ob- 
served around the other nine isolated millisecond pulsars 
in the Galactic disk. In this paper, we offer an alterna- 
tive model related to the Quark Nova (QN) phenomenon 
(Ouyed et al. 2002; Ouyed et al. 2003a). In the QN explo- 
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sion a protoplanetary disk forms from the fall-back matter. 
The disk matter gains angular momentum from the central 
quark star (QS) via the propeller mechanism. The effect 
of viscosity on the disk evolution later allow for planet 
formation. We describe these features and show how the 
points mentioned earlier are dealt with. 

This letter is presented as follows: In Sect. 2. we re- 
view the concept of QN. Sect. 3. deals with the angular 
momentum transport and the disk formation. The evolu- 
tion of the disk is described in Sect. 4 where planetary 
formation is briefly discussed. We conclude by applying 
our model to the PSR1257+12 system in Sect. 5. 

2. Quark Nova 

In the QN scenario (see Ouyed et al. 2002 and Ouyed et 
al. 2003a) the quark core of a neutron star (NS) shrinks to 
the corresponding stable compact/quark object before the 
contamination of the entire star. By contracting, and phys- 
ically separating from the overlaying material (hadronic 
envelope which is mostly made of crust material), the core 
drives the collapse (free-fall) of the left-out matter leading 
to both gravitational energy and phase transition energy 
release as high as I0 53 ergs. The QN-ejecta consist mainly 
of the NS metal-rich material outer layers. The accretion 
rate of the corresponding fall-back material is (see Ouyed 
et al. 2003a) 
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where pe is the density of fall-back matter (10 6 g/cm 3 , 
representing the neutron star crust material), R is the 
radius of the QS and M its mass. 

The rotational period of the newly formed QS can eas- 
ily spin up to millisecond periods (here we adopt a period 
of P = 2 ms) and acquire magnetic field of the order of 
B ~ 10 13 - 10 14 G. The newly born QS is defined by three 
characteristic radii: the Kcplcrian "co-rotation radius" 
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the magnetospheric radius at which the ram pressure of 
the in-falling matter balances the magnetic pressure 
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(see, e.g. Frank ct al. 1992), and the light cylinder 
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Given our fiducial values, the QS is born in the propeller 
regime, i.e. R c < i? m < Ri c (Schwartzman 1970; Illarionov 
& Sunyaev 1975), where the infalling material may be ac- 
celerated in a wind that carries away angular momentum 
from the magnetosphere and hence from the QS itself. It 



is plausible that a reasonable fraction of QNe remnants 
undergo a propeller phase after the explosion. 

The angular momentum propelled away per unit time 

is 
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The total angular momentum released in approximately 
100 seconds (the propeller lifespan; see Ouyed et al. 2003a) 
is thus of the order of 1.9 x 10 47 crg • s carried away by 
10 30 g of fall-back material. We note that a QS period can 
be as small as 0.5 ms with the limit set by the Kepler fre- 
quency (e.g., Glendenning 1997). In most cases the fastest 
spin period is smaller than 1 ms; the limit on how fast 
a gravitationally bound star (like NS) can rotate before 
being ripped off by centrifugal forces. That is, the angular 
momentum extracted by the propelled matter can eas- 
ily be as high as <~ 10 48 erg • s if P; <~ 0.5 ms. It is still 
only 5 — 10% of the total initial angular momentum of the 
QS. Almost maximal angular momentum transfer might 
be possible, but that would require specific conditions: 
very large accretion rates (and therefore a large amount 
of fall-back matter) together with a short propeller time 
to be able to compete with energy losses by gravity waves. 

Using eq.(ll) in Ouyed et al. (2003a), we estimate that 
in 100 seconds, the QS with Pi = 2 ms would have spun 
down (due to the propeller) by 30% to reach a Pf = 2.6 ms 
period. Given the initial rotational energy of the QS 

t-, „ , ,,^1 / M N , R N 2/2ms N 2 

^o,= 7.4X10 « erg (_)(—) (— ) , (6) 

this implies 2 x 10 51 erg in rotational energy is lost to the 
propelled material. 



3. Torus formation 

In order to keep some material bound to the star the 
propelled matter must first interact with the dense in- 
falling matter, releasing energy and angular momentum 
to it. Indeed, the potential energy of the propelled matter, 
—GMm pTop ./R <~ — 10 50 erg, is not small enough in com- 
parison to its rotational energy. The most realistic view 
of the process is the following: the matter propelled away 
will be subject to the infalling matter. From the numbers 
above it is straightforward to see that propelled matter 
should interact with 10-20 times more matter (from the 
surroundings) to remain bound. This mixed matter (pro- 
pelled and falling-back sharing the angular momentum) 
through viscous effects evolves into a torus. The more ac- 
curate study of the angular momentum transfer between 
the propelled matter and the surroundings would certainly 
require the use of numerical simulations and is beyond the 
scope of this paper. Here, and in order to carry on with 
our investigation, we again suggest that the outcome of 
this process is a torus or a protoplanetary disk. 
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The evolution of the torus is governed by the angular 
momentum transfer due to turbulent viscosity (Shakura 
and Sunyaev 1973) 

r 2 

v = ac s H = a -± = aVLH 2 , (7) 

where the parameter a ~ (H t / H)(v t /c s ) < 1 includes 
all uncertainties of the turbulence. Here, H t is the largest 
eddy size while H is the vertical height of the disk. The ra- 
dial diffusion timescale is r„(r) ~ r 2 /v, and if the viscosity 
is expressed by the above a formula, the viscous evolution 
time of this torus is r„(r) = a _1 (r/.ff) 2 fi _1 where f2(r) 
is the angular velocity at radius r. That is, the torus will 
expand within some seconds to a distance of 10 3 — 10 4 km 
since r/H <~ 1 (the ratio of the disk radius and the scale 
height) and assuming a ~ 0.01 — 1. 

In these early phases the accretion to the QS is large, so 
the viscous time scale for most of the torus to be accreted 
is at best some Keplerian orbital times. Nevertheless, we 
expect the star to remain bare for two reasons: i) the pro- 
peller is still efficient at deflecting matter and ii) the star's 
surface temperature is very high (> 10 McV). Therefore, 
it is natural to expect that any normal matter that man- 
aged to evade the propeller is ejected away as a hot wind 1 . 
The torus is also very dense and would likely survive the 
hot wind and any radiation field. 

4. Disk Formation and Evolution 

The torus described above will further expand to form a 
disk. The formation and evolution of such a disk is also 
governed by the angular momentum transfer due to tur- 
bulent viscosity. While most of the matter is accreted to 
the QS, most of the angular momentum is carried out with 
the part of matter demigrated to some larger distances. 

If one approximates H ~ (0.01 — 0.1)r and a <~ 0.01, 
the viscous timescale for the disk reaching r = 1 AU is 
T v( r ) ~ 10 3 — 10 5 yrs. The volume of a disk with r ~ 
10 13 cm and with a scale height H ~ 10 12 cm is roughly 
irr 2 H ~ 10 38 cm 3 . With the total mass of 10 28 - 10 29 g, 
the average density is 10" 10 — 10~ 9 g/cm 3 , corresponding 
to an average proton density of 10 14 — 10 15 cm~ 3 . The 
number density of particles is expected to be lower since 
the disk is rich in heavy nuclei. 

At later stages, the QS surface temperature would have 
decreased allowing for crust formation (channelled along 
the field line). The kinetic energy of the subsequently ac- 
creted gas may be transformed into emission of the QS 
atmosphere in the following way (e.g., Xu 2002): The mag- 
netic field channels the gas motion along the field lines to 
the QS. In this case, the kinetic energy of ions at the sur- 
face is about 100 MeV/nucleon, which is 5 times greater 
than the Coulomb barrier at the quark surface (Alcock et 
al. 1986). Accreted particles penetrate through the quark 
surface, and they are dissolved into quark matter. As a re- 
sult, the quark core is heated at the magnetic poles. The 

1 The physics for calculations of gas outflow from the crust 
is similar for both NSs and QSs (Usov 1997). 



process of heat transport through the core is very fast be- 
cause of the very high heat conductivity of quark matter, 
and therefore the quark core is nearly isothermal. Then 
the energy that is released in the process of gas accre- 
tion is radiated from the normal-matter atmosphere (the 
crust) more or less isotropically. However, this process oc- 
curs much later in the evolution, the disk then being ge- 
ometrically thin enough to survive photodissociation as 
discussed below. 



The luminosity can be estimated conservatively as fol- 
lows: if the whole disk matter is accreted to the central 
object within reversed viscous timescale, the rate should 
be at most ~ 10 29 g/10 5 yr ~ 3x 10 16 g/s, which is roughly 
0.03 of the maximal rate corresponding to the Eddington 
luminosity. When accreted to the surface of the QS, the 
total released energy in radiation is L = rjrhc 2 , where r\ is 
the accretion efficiency; in our case conservatively r\ ~ 0.15 
(it seems appropriate as a first approximation to take the 
accretion efficiency of QSs to be close to that of black holes 
and of neutron stars). The QS luminosity, the disk geom- 
etry and viscous dissipation will define the temperature 
of the disk. The dissipation can be neglected considering 
the above mentioned accretion rates (see e.g. Ruden 1993). 
The black body temperature of a perfectly absorbing plane 
with its surface inclined at an angle tan (3 = H/2r = 0.005 
to the radiation is T = [rjrhc 2 tan/3/(cr47rr 2 )] 1 ' /4 - 1330 K 
at r = 0.2 AU and less than 600 K at r = 1 AU (using 
i] = 0.15); assuming Le one obtains 3500 K and 1600 K, 
respectively. Lower temperatures are possible since for a 
rapidly rotating compact object where the accreted mat- 
ter closely corotates with the surface (the matter "softly" 
lands onto the surface), as in our case, the efficiency can 
be smaller (rj ~ 0.05; Sibgatullin & Sunyaev 2000). 

Since the baryon density of the disk is n <~ 10 14 — 
10 15 cm~ 3 , the mean free path of photons is A = 
l/(n<7x) ~ 10 9 — 10 10 cm. Therefore radiation docs not 
penetrate the inner parts of the optically thick disk, and 
it can cool to the temperatures estimated above. Added 
to that, the high energy wind particles cannot penetrate 
it either, since the Coulomb cross section is larger than 
the Thomson cross section. Therefore the disk can easily 
protect the forming dust grains (and later planetcsimals) 
from the possible radiation of the QS. Iron and most of the 
other chemical components will condense into dust grains 
when the disk has cooled to temperatures below 1500 K 
(see e.g. Lewis 1995). The condensation timescale for dif- 
ferent dust grains (e.g. Sedlmayr 1993) are roughly some 
months up to some tens of years. The dense metal rich en- 
vironment with viscous timescale that is longer than con- 
densation timescale makes the QN disk suitable for dust 
grain formation. The growth of these grains and later plan- 
etcsimals (see e.g. Lissauer 1993; Ruden 1999) will lead to 
the plausible formation of planetary bodies up to some 
earth masses, with small (up to some AU) circular orbits. 
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5. The PSR1257+12 system 

Before applying our model to the case of PSR 1257+12, 
we note that there has been some observational evidence 
of a planet or a companion orbiting the pulsar system 
PSR 1620-26, with a mass around 0.01 M@/ sin i and the 
semimajor axis of around 40 AU (Joshi & Rasio, 1997). 
This is however a triple system in a globular cluster and 
the formation history might be very different from the 
one we are considering here (see e.g. models of Sigurdsson 
1993; Ford et al. 2000). 

The three confirmed planets of the PSR 1257+12 have 
a total mass of ~ 4 x 10 28 g. Most of the disk material 
(more than 90%) in our model would have fallen back to 
the QS leaving only (0.01 - 0.1) x 10 30 = 10 28 - 10 29 g 
orbiting the QS. The total angular momentum of the 
PSR 1257+12 planets is - 1.3 x 10 48 erg • s, which is an or- 
der of magnitude higher than our estimates (see Sect. 2.). 
This can be accounted for (without stretching the model) 
by considering a longer propeller life span (>■ 100 s) 
and/or a faster spinning QS (Pi ~ 0.5 ms), which could 
then provide as much angular momentum as there is in 
the PSR 1257+12 planets 2 . 

The difficult challenge, however, is explaining 
PSR 1257+12 measured dipole magnetic field 3 which is 
B m 8.8 x 10 8 G. With such a weak field the propeller 
would not have worked making it difficult to form a 
protoplanetary disk. In our model, the importance of 
the propeller is two-fold: i) it supplies enough angular 
momentum for disk formation and ii) delays the for- 
mation of the crust thus drastically reducing accretion 
luminosity and stellar wind during the early stages of 
disk formation. One possibility is that the magnetic field 
has decayed. While this hypothesis would have difficulty 
if PSR 1257+12 is a NS (see Camilo et al. 1994 for a 
discussion) this is not necessarily the case if one assumes 
that PSR 1257+12 is a QS: an interesting feature of 
such an object is the plausible decay of the magnetic 
field due to the Meissner effect at its surface (Ouyed et 
al. 2003b), a notion which remains to be confirmed. We 
thus argue that PSR 1257+12 might be a QN compact 
remnant (born as a millisecond QS) which experienced 
such a phenomenon. Finally, comparing the QN rate of 

1 per million years per galaxy to that of the SN rate 
of 2 x 10 4 per million years per galaxy, on average 1 in 

2 x 10 4 compact objects are QN remnants. We suggest 
that a good portion of this QS population may harbour 
a planetary system as described in this work. 
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standard magnetic dipole radiation spin-down models. 
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